Hermitian four-well potential as a realization of a PT-symmetric system 
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A "PT-symmetric Bose-Einstein condensate can be theoretically described using a complex optical 
potential, however, the experimental realization of such an optical potential describing the coherent 
in- and outcoupling of particles is a nontrivial task. We propose an experiment for a quantum 
mechanical realization of a PT-symmetric system, where the 'PT-symmetric currents of a two- 
well system are implemented by coupling two additional wells to the system, which act as particle 
reservoirs. In terms of a simple four-mode model we derive conditions under which the two middle 
wells of the Hermitian four-well system behave exactly as the two wells of the PT-symmetric system. 
We apply these conditions to calculate stationary solutions and oscillatory dynamics. By means of 
frozen Gaussian wave packets we relate the Gross-Pitaevskii equation to the four-mode model and 
give parameters required for the external potential, which provides approximate conditions for a 
realistic experimental setup. 
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In quantum mechanics an observable is described by 
an Hermitian operator. This is true in particular for 
the energy, which is represented by the Hamiltonian. 
The Hermicity is sufficient for purely real eigenvalues, 
but is this really a necessary condition? Bender and 
Boettcher found that for non-Hermitian Hamiltonians 
with a weaker condition, namely VT symmetry, there 
exist parameters for which the energy eigenvalue spec- 
trum is purely real [1 , where VT stands for a combined 
action of parity V {x ^ —x, p — >■ —p), and time reversal 
T {x ^ X, p ^ —p with T'l = — iT). 

Due to the close analogy between the Schrodinger 
equation and the equations describing the propagation 
of light in structured wave guides, a PT-symmetric op- 
tical system could be visualized [2^ and experimentally 
investigated 121 H] • The necessary complex potential cor- 
responds to a complex refractive index, which is real- 
ized by balanced gain and loss of light in the wave guide. 
Several other systems with VT symmetry have been sug- 
gested and partially realized, including lasers [SHZ] , elec- 
tronics P-flD] , microwave cavities [TTl , and quantum field 
theories [U [12]. But up to date, a quantum mechanical 
realization of a PT-symmetric system is still missing. 

It was proposed [2] that a system similar to complex 
refractive indices in wave guides could be realized with 
Bose-Einstein condensates (BECs) in double-well poten- 
tials, where particles are injected in one well and removed 
from the other one. BECs in PT-symmetric double-well 
potentials have been investigated in the Bose-Hubbard 
model and the mean-field approximation |13H17| . In 
these Refs., the VT symmetry is given by a complex po- 
tential which fulfills the condition V*{x) = V{—x) and 
describes the coherent in- and outcoupling of atoms into 
and from the system. It has already been shown that a 
bidirectional coupling between two BECs is possible and 
at the same time particles may be continuously ejected 
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Figure 1. (a) PT-symmetric two- well system with tunneling 
rate J and complex potential ±ir, which describes coherent 
in- and outcoupling of atoms, (b) Alternatively, two addi- 
tional wells serve as particle reservoirs and induce particle 
currents that - under application of specific system parame- 
ters to be determined - are equivalent to the PT-symmetric 
currents of system (a). 



[iHiiTn]. 

In this letter we follow a different approach. Instead 
of injecting and removing particles, we couple two addi- 
tional wells to a double-well system. The tunneling of 
the outer wells can be used to add and remove parti- 
cles from the inner wells (see Fig. [T|). We investigate, if 
just considering the two middle wells, they behave like 
the two wells of the T'T-symmetric two-well system. For 
the four-well potential we assume a combination of four 
Gaussian beams. 



V;xt(^") = ^l^iBXp 



2x2 
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, (1) 



with Vi < the depth and si the displacement along the 
z axis of each well, and the width of a single well in 
each direction. The dynamics of a BEC is described by 
the Gross-Pitaevskii equation (GPE) 



itidtip{r,t) = 



2m 



Kxt(r)+iV5lV'(r-,i)r 



V'(r,0, 
(2) 



where g = ATrti^as/m is the strength of the nonlinearity, 
with as being the scattering length. 
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It was shown that the "PT-symmetric two-mode model 
shows the features specific of T'T-symmetric systems . 
Therefore we use the simple, but instructive two- and 
four-mode models for our investigations. For the mo- 
ment we neglect particle interaction, which will be taken 
into account later. The Hamiltonian of the linear, VT- 
symmetric two-mode model is given by 



ir -J 
-J -ir 



(3) 



which models a double-well system of non-interacting 
bosons. The real quantity J > designates the tun- 
neling amplitude between the two wells, and T G M gives 
the strength of the imaginary PT-symmetric potential, 
which models gain and loss, respectively. Obviously, this 
Hamiltonian commutes with the combined action of par- 
ity and time reversal, \PT,H^'^^] — 0, where the parity 
operator is given in this matrix representation by 



V 



1 

1 



(4) 



and time reversal is simply the complex conjugation. 

The eigenvalues and eigenvectors are easily found to be 

E± = ±V J2 - r2 and V± = (ir ± VJ^ - r2, -jf (un- 
normalized). For T < J the eigenvalues are purely real 
and the eigenvectors obey VT symmetry. In the oppo- 
site case, r > J, the VT symmetry is broken and the 
eigenvalues become purely imaginary. Thus, the simple 
two-mode model features the basic properties of a PT- 
symmetric Hamiltonian. 

To gain deeper insight, and for an easier comparison 
with the later defined four-mode model, we give the time 
derivatives of the observables of the system. Two ob- 
servables are the number of particles in each well, which 
are given by nt = i/jlipk, where G describes a 
quantum mechanical state. With the particle current 
ji2 — (V'i'02 ~ "01 V'2) between the two wells the time- 
scaled {t ht) Schrodinger equation can be brought to 
the closed set of differential equations for the observables 



dtUi 



-ji2 + 2rni, dtn2 = +ji2 - 2Tn2, (5a) 
dtJi2 = 2J^{n,-n2). (5b) 



The imaginary potential induces particle currents from 
and to the environment, jgi — 2Tni and j2e = 2rn2i 
both of which are proportional to the number of particles 
in the corresponding well. 

It is now our main purpose to investigate, whether the 
behavior of the PT-symmetric two-mode model (|3| can 
be described by a Hermitian four-mode model, where two 
additional wells are coupled to the system. The Hamil- 
tonian of the four-mode model is given by 




(6) 



The two middle wells will be symmetric, hence Ei = 
E2 — 0. The tunneling amplitudes and on-site energies 
of the outer wells are Jqi, J23, and Eq, E^, respectively. 
They may be time-dependent, as denoted by the explicit 
time-dependence in Eq. (|6|. To be able to compare the 
four-mode model with the two-mode model, we calculate 
the time derivatives of the particle populations, which 
yields the simple relations 



dtnn = -joi, 
dtn2 = J12 - ^23, 



dtni = joi - ji2, 

dtns = ^23, (7) 



where 



jki = iJki {ipk'ipt - 4'1'ipi) 



(8) 



is the particle current between adjacent wells. By addi- 
tionally considering the time derivative of ji2 , we obtain 



dtni = joi - ji2, dtn2 = ji2 - 323, (9a) 
dth2 - 2Jl^ {712 - «i) + J12 (^23^13 - J01C02) , (9b) 

where we defined Cki — ^pk4'i+^k'^i- Comparing Eqs. ([s]) 
and (|9| we can conclude that if the conditions 

joi = 2rni, J23 = 2Tn2, J01C02 = ^23^13 (10) 

are fulfilled, the two middle wells of the Hermitian four- 
mode model ^ behave exactly as the two wells of the 
PT-symmetric two-mode system ([s]). 

We now need to give the explicit time-dependency of 
the free parameters of the Hamiltonian ([6|, Eq, E3, Jqi 
and J23, such that Eqs. (10) are fulfilled. For the tunnel- 
ing elements we can set 



Jo 



dC 



13, 



J2: 



dCo 



(11) 



where d is a time-independent parameter, and can be 
tuned to bring the tunneling elements into an experimen- 
tally realizable range. The currents joi f^nd ^23 have to 
fulfill Eqs. (10 1. However, with the tunneling elements 
determined by Eqs. (11), there are no free parameters 
left to adjust these currents. Instead, we take the time 



derivatives of Eqs. (10), which yields 



a^joi = 2rni + 2rni 2rni + 2r (2rni - jia) , (12a) 

dt323 = 2rn2 + 2rn2 = 2r?i2 + 2r (ji2 - 2rn2) . (12b) 

where we allow the parameter T to be explicitly time- 
dependent. For the last equalities we used Eqs. ([t]) 
and (10). Here, T is not a quantity entering the Hamil- 



tonian directly as in the two- mode model (|3|, but a free 
parameter, which determines the matrix elements of the 
Hamiltonian (|6|. 

Now we can calculate the time derivatives of joi and j23 
from the definition (|8|. This leads to the linear system 
of equations for the on-site energies Eq and E3, 



JoiCqi djoiji3 
-djo2323 —J23C23 



Eo 
E3 



(13) 
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Figure 2. (Color online) (a) Populations rii = n2 — 1/2 
of the middle wells, which are both equal and constant in 
time, (b) Populations no (solid) and na (dashed) of the outer 
wells, which depend linearly on the time, i. e. no ~ — P and 
ns = P. (c) Time-dependent on-site energies Eq (solid) and 
Es (dashed), (d) Tunneling amplitudes Joi (solid) and J23 
(dashed). Time is given in units of to — h/ J\2- 



with the entries 

vo = 2tni + 2T (joi - ji2) - 2 Jq^ (no - ^^i) 

-J01J12C02 - d ( Joijo3 + -^12^23 - J23jl2) joi, (14a) 

V3 = 2tn2 + 2T (ji2 - j23) - 2 JI3 (n2 - 713) 

+ -^23-^12Ci3 - d ( J01J12 - Jl2joi - •^23^03) ^23- (14b) 

Here we have defined the modified currents jj./ — 
i{ipktpi — tpki^i)- Thus, the on-site energies Eq and 
are used to maintain the validity of Eqs. (10 1. Since the 



energies do not determine the currents but their time 
derivatives, the initial wave function must be chosen in 
such a way that the conditions are fulfilled. 

So far we could find the explicit time-dependences of 
the matrix elements Eo{t), Ez{t), Joi(i), and J23(i) of the 
four-mode model in order that at every time the two mid- 
dle wells have the same behavior as the T'T-symmetric 
two-mode model. Thus our method is a valid possibil- 
ity to realize a T^T-symmetric quantum mechanical sys- 
tem. To calculate the matrix elements at every time step 
we integrate the four-dimensional complex Schrodinger 
equation with a numerical integrator and use Eqs. (11) 
and (13 1. We now give two examples for different solu- 



tions. 

First we consider a quasi-stationary solution, i. e. a 
state in which the particle numbers in the two middle 
wells are stationary. These states correspond to the sta- 
tionary solution of the two-mode model. We prepared 
this stationary solution for T / J12 = 0.5 at i = 0. Fig. [2] 
shows the results. As required, the number of particles in 
the middle wells are equal and constant in time (we have 
chosen the normalization such that ni + 712 = 1). From 
Eqs. ([7| and ( 10 ) we then obtain no — — F and — F, 
plotted in Fig^'b). Figs. [2j^c) and (d) show the calcu- 
lated matrix elements. All of them vary only slightly in 
time. Due to the linear decrease of the number of parti- 
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Figure 3. (Color online) Same as Fig.|2] but for oscillatory dy- 
namics. Additionally, the total number of particles in wells 1 
and 2 are plotted as black crosses in (a). 



cles in well 3 the available time for the VT symmetry is 
limited. 

As a second example we prepared a non-stationary so- 
lution at t = for r/Ji2 = 0.5 (Fig. |3|. There are the 
typical Rabi-type oscillations, but with a smaller phase 
difference A(/) < tt, which leads to a non-constant added 
number of particles in the two middle wells. This value 
oscillates harmonically around its mean value. The same 
behavior is obtained for the optical system in Ref. [2]. 
The matrix elements show a quasi-oscillatory behavior. 
As before, the time for an exact VT symmetry is limited. 

For these calculations, the initial conditions had to be 
chosen appropriately so that they obey VT symmetry. 
This would be a difficult task in an experiment. For 
that, we propose an approach of adiabatically increasing 
the VT parameter F. We start at the ground state of 
the four-well system, and increase F according to T{t) = 
Ff [l-cos(7rt/tf)]/2 for t G [Q,tf\ (see Fig. |4|i)). The 
quantities Ff and tf have to be chosen such that F/F <C w, 
where uj is the typical oscillation frequency. The system 
then changes adiabatically from the ground state of the 
Hermitian system to a T'T-symmetric ground state. 

So far we have neglected the contact interaction of 
the atoms. However, as we show next, taking into ac- 
count the nonlinear contact interaction will not change 
the characteristic T'T-symmetric behavior. The Hamilto- 
nian, which is nonlinear in the mean-field approximation, 
then reads 

H^'\t) = H[^{t) + cdiag [m^ , 1V2I' , , \u 

(15) 



where H^^^ is the linear part ([6|. The quantity c mea- 
sures the strength of the interaction. We have to recal- 
culate the time derivatives of the observables. For the 
derivatives of the populations, rii and 71-2, we obtain the 
same as in the linear case. For the current between the 
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middle wells we obtain 

Two wells: dj = {dtj)ii,-, + Jc{ni - nj) C12, (16a) 

Four wells: dtji2 = {dtji2)iin + J12C (ni - ^2) C12. 

(16b) 

Generally, the time evolution of C12 differs for the two- 
and four-mode model. But in the case of adiabatically 
increasing F, we have ni « 712, i.e., with the choice of 
parameters given above we have an approximate equiva- 
lence of the two- and four-mode model. Solely the linear 
system of equations ( 13 ) has to be modified to include 
the interaction. 

We further want to give approximate parameters for a 
realistic potential. Four wells can be realized by a super- 
position of four Gaussian laser beams ([ij . The dynamics 
of a BEG is described by the GPE (|2|. To relate the four- 
mode model to the GPE, we assume the wave function 
to be a superposition of frozen Gaussian wave packets, 



(17) 



k=0 



For simplicity we assume the widths to be constant in 
time, and the same for each Gaussian, and the displace- 
ment of each Gaussian to be the same as the displacement 
of the corresponding well, = . Only the amplitudes 
of the Gaussians d*^, are dynamical variables. 

Multiplying Eq. (|2| by 1/)* and integrating over R^, we 
obtain the equations of motion for d'^, 



k=0 



/c=0 



H 



(18) 



With the method of symmetric orthogonalization [H] we 
can write these equations as a Schrodinger equation with 
a symmetric and Hermitian 4x4 Hamiltonian. By con- 
sidering only nearest neighbors in the integrals, we can 
relate the matrix elements of the four-mode model ^ to 
the realistic potential ([l]). This yields 



2m 



Ay+A, 



VkPxMz, 



Jlk 



2m 



_ ^1/{1+A,w',) 



AxAyAz 



m V TT 
with the abbreviations 



/3a 



1 + A^wl ' 



7 = exp 



' 2 



(19a) 

(19b) 
(19c) 

(20) 



The parameters A^ are determined by minimizing the 
mean-field energy. 
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Figure 4. (Color online) (a) Adiabatic ramp of the current 
from r = to a target value of Ff as a function of the time t. 
Time dependences of (b) the potential depths Vb (solid) and 
V3 (dashed) and (c) the displacements of the outer wells (5o 
(solid) and ^3 (dashed), (d) and (e) Number of particles in 
each well. 



We calculated the adiabatic current ramp for a con- 
densate of = 10^ atoms of ^^Rb with a scattering 
length tuned to a = 10.9aB (ae being the Bohr radius). 
The distance between the middle wells \& I — 2 |j,m. We 
express the distance of the outer wells by their devia- 



tion from a equidistant lattice, i.e. 

-.3, 



31/2 + 6° and 
2 /^;2 



= 31/2 + 6^. The basic unit of energy is Ei = It? /ml 
which yields Ejh — 29.1 Hz. As initial condition we use 
the ground state for Vi — V2 — —8OE1 and Vq = V3 = 
— \22Ei. The parameter d in Eq. ( 11 ) is chosen such that 
5^ — 5^ = Q for t < 0. We have chosen the widths of the 
trap (jlj to be = Wy = 41 and Wz = 1/2. The basic 
unit of time is t; — h/Ei = 5.47 ms. Fig. |4] shows the 
results for T j / J12 = 0.5 and tf/ti 70. 

The trap depth V§ has to be increased, whereas V§ 
has to be decreased, in an almost symmetric way. The 
distance of the outer wells has to be varied within a few 
percent of the distance /. The results of the particle num- 
bers in each well confirm that the change of F is approx- 
imately adiabatic (see Fig. |4]). For these conditions, the 
system arrives at an approximate T^T-symmetric ground 
state for t > 70ti, which can be determined by the con- 
stant number of particles in wells 1 and 2. We note that 
the change of the distances of the outer wells can be ne- 
glected, as further calculations show. Furthermore, the 
simulation indicates that the system is robust with re- 
spect to small random perturbations of the external po- 
tential. 

To summarize we have shown that the two middle wells 
of the Hermitian four-mode model can show the same be- 
havior as the T^T-symmetric two-mode model, and thus 
offers an approach to realize VT symmetry in a quan- 
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turn mechanical system. This agreement is exact in the 
absence of interaction. We have proposed the method 
of adiabatically increasing the VT parameter F to cre- 
ate the T'T-symmetric ground state. This works also 
approximately for particles with interaction in the mean- 
field limit . We finally estimated parameters of a realistic 
potential, which would be necessary to prepare such an 
experiment. The time-dependent potential ([T|) could be 
realized e. g. using an acousto-optical modulator [22 . For 
future work it is desirable to extend these investigations 



to a BEC described by the full GPE to obtain more ac- 
curate parameters for an experimental realization. 



ACKNOWLEDGMENTS 

This work was supported by DFG. M. K. is grateful 
for support from the Landesgraduiertenforderung of the 
Land Baden- Wiirttemberg. 



[1] C. M. Bender and S. Boettcher JPhys. Rev. Lett.| 80, 5243 
(1998) 



[2] S. Kl aiman, U. Giinther, and N. Moiseyev, Phys. Rev. 



Lett 



_101, 080402 (2008) 
Guo, G. J. Salamo, D. 



Duchesne, R. Morandotti, 
M. Volatier-Rav at, V. Aimez, G. A . Siviloglou, and D. N. 
Christodoulides, |Phys. Rev. Lett. 1 103, 093902 (2009) 
[4] C. E. Riiter, K. G. Makris, R. El-Ganainy, D. N. 
Christodouhdes, M. Segev, and D. Kip, Nat. Phys. 6, 
192 (2010) 

[5] Y. D. Chong, L. Ge, and A. D. Stone, Phys. Rev. Lett.| 

106, 093902 (2011) 
[6] L. Ge, Y. D. Chong S. Rotter, H. E. Tiireci, and A. D. 

Stone, iPiyl Rev. A|84, 023820 (2011) 



[7] M. Liertzer, L. Ge, A. Cerjan, A. D. Stone, H. E. Tiireci, 

and S. Rotter, Phys. Rev. Lett. 108, 173901 (2012) 
[8] J. Schindler, A. Li, M. C. Zheng, F. M. EUis, and T. Kot- 

tos, Phys. Rev. A| 84, 040101 (2011) 
[9] H. Ramezani, J. Schindler, F. M. Ellis, U. Giinther, and 
T. Kottos, jPhys. Rev. A, 85, 062122 (2012) 
[10] J. Schindler, Z. Lin, J. M . Lee, H. Ramezani, F. M. Ellis, 

and T. Kottos, J. Phys. A 45, 444029 (2012) 
[11| S. Bittner, B. Dietz, U. Giinther, H. L. Harney, M. Miski- 



Oglu, A. Richter, and F. Schafer, |Phys. Rev. Lett.| 108, 

024101 (2012) 

[12] C. M. Bender, V. Branchina, and E. Messina, ' Phys. Rev. 
r ^85, 085001 (2012) 

[13] E. M. Graefe, U. Giinther, H. J. Korsch, and A. E. 

Niederle, J. Phys. A 41, 255206 (2008) 

[14| E. M. Graefe, H. J. Korsch, and A. E. Niederle, j| Phys. 
I Rev. Le tt. 101, 150408 (2008) 

jl5| E.-M. Graefe, H. J. Korsch, and A. E. Niederle, |Phys. 
[~ Rev. A 82, 013629 (2010) 

\T6\ H. Cartarius and G. Wunner, Phys. Rev. A| 86, 013612 

(2012) 



[17] H. Cartarius, D. Haag, D. Dast, and G. Wunner, J. Phys. 
I Al45, 444008 (2012) 

[18] "Y. Shin, G. B. Jo, M. Saba, T. A. Pasquini, W. Ketterle, 
and D. E. Pritchard, Phys. Rev. Lett. 95, 170402 (2005) 

[19| R. Gati, M. Albiez, J. Folhng, B. Hemmerling, and M. K. 
Oberthaler, AppUed Physics B 82, 207 (2006), ISSN 
0946-2171 

[20] E.-M. Graefe, [TPhys. A|45, 4 44015 (2012) 
[21] P.-O. L6wdin,<J. Chem. Phys.118, 365 (1950) 
[22[ K. Henderson, C. Ryu, C. MacCorinick, and M. G. 
Boshier, ,New J. Phys. 11, 043030 (2009) 



